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SUMMARY 

A  method  of  analysis  is  presented  hy  means  of  which  certain  missile 
parameters,  in  particular  the  vacuum  specific  Impulse,  may  be  oaloulated 
using  only  the  telemetry  record  of  longitudinal  acceleration.  The  accuracy 
of  the  method  is  discussed  on  a  statistical  basis  and  graphs  are  provided 
which  may  be  used  to  predict  the  accuracy  of  calculation  for  a  missile  of 
loiown  nominal  performance. 
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1  INTRODUCTION 

Dxiring  the  course  of  the  flight  trials  of  a  ballistic  missile,  it  is 
necessary  to  measure  the  performance  of  the  rocket  engine.  The  specif i- 
oation  limits  for  most  ballistic  missiles  are  so  narrow  as  to  demand  a  high 
order  of  accuracy  of  measurement  of  engine  performance.  One  of  the  most 
important  parameters  is  the  specific  impulse  and  the  usefulness  of  the 
method  presented  in  this  Note  lies  in  the  fact  that  it  allows  the  specific 
impulse  to  be  derived  directly.  A  considerable  part  of  the  Note  is  oonoemed 
with  the  accuracy  of  the  method  and  graphs  are  given  which  may  be  used  to 
predict  the  acoxiracy  of  measxirement  for  a  missile  of  known  nominal 
performance. 

2  APPARENT  ilCCELSRATION  DT  VACUO 

For  missiles  of  range  longer  than  a  hundred  miles  or  so  the  major 
part  of  the  boost  phase  takes  place  outside  the  earth's  atmosphere. 

In  this  note,  it  will  be  assumed  that  the  missile  is  in  vacuum  when 
the  aerodynamic  forces  are  negligibly  small  ownpared  with  the  thrust  of  the 
motor  so  that  neglecting  these  forces  will  produce  an  acceptably  small 
error. 


Outside  the  atmosphere  only  two  forces  act  on  the  missile,  one  being 
the  motor  thrust,  the  other  being  that  of  gravitational  attraction.  The 
missile's  true  acceleration  is  due  to  both  these  forces  acting  together 
but  to  understand  its  behavioxir  more  fxilly  it  is  useful  to  resolve  the 
'true*  acceleration  into  (i)  an  'apparent'  acceleration,  due  to  motor 
thrust  alone,  and  (ii)  a  'gravitational'  acceleration,  due  to  the  gravity 
field  alone.  In  general  these  will  aot  in  different  directions.  If  an 
accelerometer  is  mounted  in  the  missile,  then,  dxiring  its  flight  in  vacuo, 
it  will  measxire  only  the  apparent  acceleration.  This  is  because  accelero¬ 
meters  meastire  acceleration  by  sensing  the  force  which  has  to  be  exerted  on 
a  known  mass  in  order  that  the  mass  moves  with  the  same  acceleration  as  the 
missile.  However,  gravitational  attraction  exerts  forces  on  the  mass  in 
the  accelerometer  and  on  the  missile  mass  so  that  the  accelerations  due  to 
this  cause  are  the  same.  Because  there  is  no  relative  acceleration  between 
the  missile  and  the  accelerometer  arising  from  the  earth's  gravitational 
field,  no  force  oan  arise  in  the  accelerometer  mounting  system,  henoe 
a  missile-bome  accelerometer  cannot  measure  that  component  of  the  missile's 

true  acceleration  which  is  due  to  gravity.  Normally  accelerometers  have 
only  one  degree  of  freedom,  i.e.  they  measure  the  component  of  the  apparent 
acceleration  in  one  direction,  usually  fixed  in  relation  to  the  missile  body 
axes.  However,  if  the  missile  employs  inertia  navigation,  then  the  accelero¬ 
meters  will  be  moianted  on  a  table  which  will  maintain  a  fixed  direction  in 
space  independent  of  missile  orientation. 

3  EQUATION  OF  MOTION  ALONG  THE  FLI&HT  PATH 


If  f  is  the  apparent  acceleration  along  the  flight  path,  then  the 
equation  of  motion  is  at  any  instant 


f  =  Is 

m 

where  T  is  the  component  of  thrust  along  the  flight  path  (lb  wt) 
m  is  the  mass  of  the  missile  in  lb 

2 

f  is  apparent  acceleration  in  ft/seo 
g  is  52.2  ft/seo^. 


(1) 
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It  will  be  assximed  that  the  thrust  and  the  propellant  oonsumption,  m,  of  the 
engine  are  constant.  If  at  time  t  =  0,  the  mass  of  the  missile  is  m^,  then 

the  mass  of  the  missile  at  time  t  will  be  less  due  to  the  oonsimption  of  the 
propellants  and  oan  be  written 


m 


m  "“A  tj 
0  1 


therefore  f  =  - —  • 

m  -  m  t 
0 

By  virtue  of  the  above  assumption,  the  speoifio  impulse,  I,  which  is 
defined  by  the  ratio  of  thrust  to  mass  flow,  will  also  be  oonstant.  Henoe 


f  a 


JS- 


”c/“  “ 


3 


(2) 


where 


=  m/m 


We  assume  that  the  apparent  aoceleration  f  measured  aboard  the  missile 
is  telemetered  back  to  the  ground.  We  also  assume  that  the  apparent  aocelera¬ 
tion  is  sampled  along  with  other  measurements  so  that  the  output  from  -tlie 
grovmd  i«oeiver  is  a  series  of  readings  made  at  equal  intervals  of  time. 

During  the  process  of  transmission  from  the  missile  to  the  ground  reoeivers 
and  of  oonversion  to  a  suitable  form,  the  data  will  suffer  a  oertain  amount 
of  degradation.  Furthermore,  the  aocelerometer  will  respond  to  a  oertain 
extent  to  missile  vibration.  Consequently  the  data  whioh  is  available  is 
subjeot  to  error.  However,  we  shall  assume  that,  at  least  for  the  range  of 
data  over  whioh  the  analysis  is  to  be  conduoted,  the  output  from  the  ground 
reoeivers  is  a  linear  function  of  the  apparent  acceleration.  Henoe 

f  =  k(y  -  a) 

where  k  is  the  oonstant  of  proportionality 
a  is  the  zero  or  bias  error 
and  y  is  the  groxind  receiver  output. 

Substituting  for  f  in  equation  (2)  we  get 


y 


a  + 


(3) 


where 


0  =  Ig/k 


(4) 


The  method  given  below  proposes  that  the  coefficients  'a*,  b  and  o 
may  be  determined  by  fitting  a  curve  of  the  form  of  equation  (3;  to  the 
telemetry  readings,  y,  and  the  corresponding  values  of  time,  t.  Normally 
•a*  is  measured  during  the  pre-flight  calibration  of  the  aocelerometer  and 
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telemetry  ayatem.  However,  any  biaa  errora  ariaing  in  the  telemetry  eyatem 
during^  flight  will  alter  the  value  of  *a’.  It  may  he  deairahle,  therefore, 
to  determine  ’a’  hy  the  proceaa  of  curve  fitting.  Thua,  both  the  method  of 
determining  all  three  coefficient a,  ’a*,  h  and  o,  and  the  method  of  deter¬ 
mining  h  and  o,  given  the  value  of  ’a*,  to  aome  eatimated  accuracy,  are 
given  below.  The  former  oaae  is  referred  to  aa  the  caae  of  ’a*  lonknown 
and  the  latter  as  that  of  ’a’  known. 

The  value  of  c  determined  by  the  propoaed  method  may  be  uaed  to 
determine  the  specific  impxilse,  I,  by  subatitution  in  the  formula 


where  the  scale  factor,  k,  will  be  determined  to  some  estimated  accuracy 
by  the  pre-flight  calibration. 

The  coefficient,  b,  may  be  uaed  to  determine  the  mass  ratio,  X,  which 
is  defined  as  the  ratio  of  the  all-up  weight  at  launch  (t  =  O)  to  the  weight 
at  some  time  t  after  launch.  Thua 


X  = 


The  mass  at  time  t  may  be  evaluated  if  the  mass  at  laimoh  is  known. 

Thua 


m  =  m  -  mt 
o 

=  -  tA)  • 

For  a  given  accuracy  of  measurement  of  b  and  m^,  the  accuracy  of  determination 

cf  m  Td.ll  reduce  linearly  with  time,  so  that  unless  very  accurate  determina¬ 
tion  of  both  m^  and  b  are  poasibie,  the  acoxiracy  of  determination  of  m  toward 

the  end  of  the  boost  phase  is  unlikely  to  be  satisfactory. 

Alternatively,  the  mass  at  time  t  may  be  determined  in  conjunction  with 
an  independent  measurement  of  mass  flow,  m,  e.g.  by  flowmeter.  Thus 


m  = 


m  -  mt 
0 


=  fi(b  -  t) 


The  accuracy  of  determining  the  mass  by  this  method  will  not  be  dependent  on 
time  so  that  althoiagh  the  accuracy  of  measiarement  of  m  in  flight  may  not  be 
aa  good  as  the  accuracy  of  measurement  of  m^  just  before  launch,  the  latter 

method  is  likely  to  be  more  satisfactoiy  than  the  former  in  the  later  stages 
of  boosted  flight. 
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An  Independent  measurement  of  the  mass  flow,  m,  is  also  required  to 
determine  the  thrust  in  conjunction  with  the  measured  value  of  the  speoifio 
impulse,  using  the  formula 

T  =  I  m  . 


This  method  of  deriving  the  thrust  serves  as  an  independent  and  overall 
cheok  on  the  engine  performance  becavise  it  requires  no  prior  knowledge  of 
the  engine  performance  oharaoteristios. 

k-  THE  METHOD  Qg  DETERMININ&  THE  COEPFICIENTS 


The  first  oase  to  he  dealt  with  is  the  general  oase  where  ’a'  is 
\mknown.  The  more  partiovilar  case  where  'a'  is  known  may  then  be  obtained 
along  similar  lines  (see  seotion  7)« 

It  will  be  assumed  that  although  y  will  be  siibjeot  to  errors,  measure¬ 
ment  of  time  will  be  sufficiently  aoo\irate,  and  as  a  general  rule  this  is  so* 
In  general  a  large  number  of  observations  of  y  and  t  will  be  made  so  that  the 
problem  is  to  determine  the  values  of  the  ooeffioients  a,  b  and  o  whioh  best 
fit  the  observations.  This  will  be  done  by  the  method  of  least  squares' • 
Because  the  function  that  must  be  fitted  to  the  data  is  non-linear  in  t,  it 
is  necessary  to  use  a  process  of  differential  adjustment  in  order  to  deter¬ 
mine  the  ooeffioients.  The  method  has  therefore  been  treated  in  detail. 

At  a  time  t^^,  let  the  fitted  value  of  y^  be  Thus  when  a,  b  and  o 

have  their  most  probable  values,  is  given  by 


a  + 


b  -  t^ 


The  residual,  v^^,  at  time,  t^^,  is  defined  as  the  difference  between 
the  observed  value  y^^  and  the  fitted  valixe  Yj^,  i.e. 


Vi  =  yi  -  \ 


therefore 


(4) 


The  method  of  least  squares  shows  that  if  the  true  errors  (i.e.  the 
differences  between  the  true  and  the  observed  values  of  y)  are  iiidependent 
and  nozTnally  distributed,  the  most  probable  values  of  a,  b  and  o  are  those 
whioh  make  the  sum  of  the  squares  of  the  residuals  a  miMmiim,  i.e. 


1 
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These  three  equations,  which  are  called  the  normal  (or  normalised) 
equations,  are  sufficient  to  determine  a,  h  and  o.  Thus  after  squaring  (4) 
and  summing  over  all  values  of  i  and  then  partially  differentiating  with 
respect  to  a,  h  and  c  in  turn  we  get 


The  values  of  the  coeffioients  which  satisfy  the  above  equations  are 
the  most  prohahle  values.  However,  it  will  be  seen  that  the  equations  are 
non-linear  in  h  and  can  only  he  solved  hy  a  process  of  successive  approxi¬ 
mation.  Because  it  is  desirable  to  oaloulate  the  variances  of  the  coeffi¬ 
oients,  expressions  for  whioh  have  only  been  derived  for  linear  normal 
equations,  it  is  necessary  to  linearise  equation  (4)  by  the  process  of 
differential  adjustment  rather  than  to  solve  equations  (6)  direot!ly. 

The  process  of  differential  adjustment  assumes  that  good  estimates  of 
b  and  o,  say  b^  and  c^  respectively,  can  be  obtained.  Next  the  most  probable 

values  of  adjustments  3  and  y,  which  must  be  made  to  b^  and  c^  can  be  found, 

i.e.  the  values  of  a,  p  and  y  can  be  fo\and  whioh  make  the  svim  of  the  squares 
of  the  residuals  a  minimum.  It  is  supposed  that  the  original  estimates  b^ 

and  o^  are  sufficiently  good  that  p  and  y  oan  be  regarded  as  small  quantities. 
In  equation  (4)  putting  b  =  b^  -  P  and  0  =  o^  +  y  we  get 


If  and  Tg — are  small  compared  with  1,  the  above  equation  becomes 

o  °o  ^i 
approximately 

^i 


Putting 
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Squaring  equation  (7)>  summing  for  all  available  values  of  i  and 
partially  differentiating  with  sTespeot  to  a,  p  and  y  in  turn  we  get 


1 

n 


n  2  n  n 

0,  i  c 


,  'i  T  ^  1 


n 


Et‘'ES-Ei-E('.'®t-» 


(8) 


(9) 


(10) 


The  solution  of  these  equations  in  determinantal  form  is 


ijli 

and  A  is  the  oo-faotor  of  the  element  in  the  p 
PG 

determinant  D. 


row  and 


oolumn  of 


the 


Having  determined  the  values  of  a,  p  and  y,  p  and  y  oan  be  used  to 
adjust  bg  and  0^  and  the  adjusted  values  used  as  a  new  estimate  of  b  and  0. 
Pvorfcher  adjustments  oan  then  be  computed  by  the  above  method  until  b  and  0 
have  been  oomputed  to  the  desired  accuracy. 
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It-  is  not  proposed  to  consider,  analytically,  the  convergence  of  the 
method.  However  on  all  the  occasions  when  the  method  has  been  used,  the 
convergence  has  been  found  to  be  satisfactory. 


5  ESTIMATES  OF  THE  VARmCE  OF  THE  COBFFICIEM’S 

Formulae  for  the  variance  of  the  coefficients  have  been  derived^  on 
the  assumption  tliat  the  errors  are  independent  and  normally  distributed  and 
that  the  distribution  is  independent  of  time.  In  the  notation  of  the 
previous  section  these  formulae  are:- 


Tfhere 


2 

s 


2 

s 

a 


2 


2 


2 

s 


2  ^2 

®  D 
D 


1  E 
m  -  3  H 


n 


1 


m  -  3 


(15) 

(16) 
(17) 


(18) 


where  s„ 
a 

2 

2 

m 


estimate  of  the  variance  of 
estimate  of  the  variance  of 
estimate  of  the  variance  of 
niunber  of  observations 


a 

b 

0 


and 


(19) 
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It  shovild  be  noted  that  the  above  estimates  of  the  variance  relate 
strictly  to  P  or  y  rather  than  to  b  or  o.  However  as  b  =  b^  -  P  and  b^  is 

a  constant,  then  the  variance  of  b  is  the  same  as  the  variance  of  P* 
Similarly  the  variance  of  o  is  equal  to  the  variance  of  f 

6  .  PREDICTION  OF  THE  VARIANCE  OP  THE  COEFFICIENTS 


It  will  be  seen  from  eq^tion  (lO  that  the  determinant  D  and  conse¬ 
quently  the  oo-faotors  depend  only  on  o^,  b^  and  t^  and  is  independent  of  y. 

It  is  therefore  possible  to  e:q>ress  the  variance  of  the  coefficients  in 
terms  of  these  variables  after  maldLng  certain  simplifying  assumptions.  It 
is  then  possible  to  predict  the  variance  of  the  coefficients  approximately 
for  a  given  missile  and  so  obtain  prior  knowledge  of  the  precision  of  the 
analysis . 


In  the  determinant  D 


there  are  only  5  terms,  i»e» 


Z 


E 


1 

It  Should  be  noted  that  o  is 

0 


n 


1 


a  common 


factor  in  the  second  column  and  so  may  be  brought  outside  the  determinant* 
Also,  by  interchanging  the  last  two  columns  and  then  the  last  two  rows,  the 
determinant  becomes  symmetrical  about  the  leading  diagonal. 


We  now  assume  that  the  observations  are  made  at  uniform  intervals  of 
time  6t  and  that  the  Interval  is  small.  Consequently 

t 

n  n  n 


dt 

(i  - 


n 


'  *0^ 


for  p  =  0 
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Putting  =  f\i  we  get 

n 

E’ 


n 


X 


■fe  ’I 


h- 1] 


26t  X 

0 

— ^  [n^  -  1]  . 

36t 


Slibstituting  for  the  above  in  the  determinant  D  we  get 


D  = 


X  (ti  -  1) 
n^  '  ' 

6t 


o^(ti  -  1 ) 


OqCt)  -  1) 


STT—  -St 


-  1)  0  (n^  -  1) 


btx” 


6t  33{^ 


6t  2x^ 


"St 


OqC'H  “  0  (ti  -  i) 


6t  2x: 


Tt 


Expanding  this  determinant  in  terms  of  the  elements  of  the  first  row  and 
their  oorresponding  co-faotors  we  get 

.  2 


D  = 


6t' 


0  0 

^  log  ^  pLiilKsL-iLJl  .  laL^ 


2=4 


3x? 


log  n 


']] 


(20) 
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Similarly  the  oo-f actors  of  the  diagonal  elements  are 


Sl  6. 2  ■  V 


-  1)^- 


t  sT  L 
0 

o2  /„ 


Lu-iJl 


6t2  12 

0 


(21) 


(22) 


^ 


,r).(.T^.  --il .  - 1^12 


3ti 


(ri-  D' 


(23) 


Now 


1  =  :ii-i  .1 
^  ^0 


X  -  X 
n  0 


+  1 


therefore 


1  -  n  _  n  0 

•n 


therefore 


X  =  (t  -  t  )  ' 

0  '  n  o'  r\  -  ^ 


(24) 


Substituting  for  x^  in  equations  (21 )  to  (23)  and  then  substituting  for 
A,i,  ^22  and  in  equations  (I5)  to  (17)  we  get 


W-T-) 

^  n  0' 


jXztllf 

V  12ti 


(25) 


(n-l) 


ijtU; 


[log  riT 


6t(Vto)' 


(26) 


and 


where  C 


"0  '■ 


1  _  'ij, 

6t(t  -t  }  ■* 


li  (rP-0  -  T^-n-l) 


s  J  6t(t^-to) 


(27) 


”  (■n”0(Tp-l)  logTi  -  [log  T^]2j 
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It  vd.ll  be  seen  that  the  right  hand  sides  of  the  above  three  equations 
are  functions  of  7}  only.  The  left-hand  sides  contain  the  standard  deviations 
of  the  coefficients  s^,  s^  and  s^  together  with  quantities  which  are  known 

or  which  can  be  arbitrarily  chosen.  Thtzs  the  parameters  on  the  left-hand 
side  may  be  conveniently  plotted  as  functions  of  t]  as  in  Figs.1  to  3* 

In  predicting  the  standard  deviation  of  the  coefficients,  7}  may  be 
interpreted  as  the  local  mass  ratio  defined  by 


mass  of  missile  at  start  of  the  analysis 
^  “  mass  of  missile  at  end  of  the  analysis 

m  -  m*.  b  -  t  X 

therefore  n  =  — ^ ^  =  r - =  — 2  ,  (28) 

•  D  —  X  X 

“o  “  %  n  n 

It  is  assumed  here  that  the  nominal  values  of  m  ,  fi  and  o  will  be 

0  0 

used  when  predicting  the  accuracy  which  may  be  obtained  so  that  the  assump¬ 
tion  of  constant  mass  flow  throughout  the  whole  of  the  boosted  flight  will 
legitimately  apply. 

2 

The  quantity  s  vd.ll  depend  on  errors  from  all  sources,  e.g.  reading 
error,  radio  propagation  errors,  transducer  errors  eto.  However,  prior 
experience  with  similar  telemetry  and  transducer  installations  will  afford 
information  on  the  likely  value  of  s  which  may  be  realised. 

7  CASE  TOEBE  «a«  IS  KNOM 


In  this  case  and  in  the  notation  of  the  previous  section,  we  may  write 


y±  =  ^1  -  ^ 


and 


Y’  =  Y.  -  a 

1  1 


Because  'a'  is  known,  y^  will  be  known,  so  that  equation  (4)  becomes 


(29) 


Putting  b  =  b^  -  p  and  c  =  0^  +  y,  the  above  equation  reduces  to 


V 


1 


(30) 


Squaring  this  equation  and  sviinming  for  all  available  values  of  i,  and  then 
partially  differentiating  with  respect  to  P  and  y  in  turn  we  get 
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n  n 


where,  as  before,  x.  =  b  -  t . . 

*  ’101 

The  solution  of  these  equations  is 


(31) 


(32) 


where 


The  method  of  repeated  adjustment  of  b  and  0  is  the  same  as  in  the 
general  case .  ° 

The  estimate  of  the  variances  of  the  coefficients  are  given  by 


n 
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n 


(39) 


Making  use  of  the  approximations  derived  in  section  6  for  the  summations, 
the  expressions  which  may  he  \ised  to  predict  the  accuracy  of  the  coefficients 
are  in  this  case 


% 


M 


12t1^ 


(W) 


1°  I  ,  ^ 

s  J  5t(t^  -  t^) 


(M) 


The  left-hand  sides  are  the  same  as  the  parameters  on  the  left-hand 
sides  of  the  corresponding  equations  (26)  and  (27)  in  the  general  case. 
Similarly  the  right-hand  sides  are  functions  of  t).  Equations  (W))  and  (41 ) 
are  represented  graphically  in  Figs. 4  and  5* 

The  variances  of  h  and  0  given  by  equations  (40)  and  (41)  apply  only 
when  'a*  is  known  precisely.  In  practice,  the  pre-flight  measurement  of  'a* 
will  be  subject  to  errors,  while  further  bias  eiror  effectively  altering  'a' 
will  occur  in  flight.  Then  error  in  the  knowledge  of  'a'  will  propagate 
errors  in  both  b  and  c.  Let  the  error  be  6a.  Therefore  yj^  =  y^  -  a  -  6a. 

Thus  the  error  propagated  in  b  and  c  may  be  fo\md  by  substituting  this 
expression  for  yj^  in  equations  (53)  and  (34)  and  separating  out  that  part 

containing  6a,  Thus 
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D» 


Making  use  of  the  approximations  derived  in  section  6  these  equations 
reduce  to 


&a 


(t  - 1  r 

^  n  o' 


=  r(n-i)^  -  K'n^-O  log  Ti 


(n-1)' 


'J 


(Mf) 


6r  _  n 

"  (-n-1)^  I 


log  Ti  -  6(n-l)  .  {hd) 


These  equations  are  represented  graphically  in  Pigs. 6  and  7, 


8  DISCUSSION 


In  practice,  the  usefulness  of  the  method  of  determining  the  missile 
parameters,  b  and  c,  will  depend  mainly  on  the  accuracy  with  which  they 
can  be  found.  This  section,  therefore,  will  be  concerned  chiefly  with  an 
appreciation  of  those  factors  which  have  a  significant  effect  on  the 
aoouraoy. 

The  curves  given  in  Pigs.1  to  7  may  be  used  to  assess,  approximately, 
the  aoouraoy  which  can  be  obtained  by  applying  the  method  of  analysis  to  a 
missile  of  known  nominal  performance.  Although  these  curves  are  of  general 
application,  they  do  not  show  clearly  how  the  aoouraoy  is  affected  by  the 
time  of  the  start  of  the  analysis,  t  ,  and  the  period  of  the  analysis 

t^  -  t^.  This  arises  from  the  faot  that  the  local  mass  ratio,  ti,  is  depen¬ 
dent  upon  both  t^  and  t^  -  t^.  Consequently,  the  sensitivity  to  t^  -  t^  is 

hidden,  to  a  certain  extent,  by  the  choice  of  n  as  the  independent  variable 
in  Figs.1  to  7*  In  order  to  illustrate  this  sensitivity.  Pigs. 8  and  9  have 
been  prepared  as  a  particular  example.  These  figures  show  the  variation  of 
the  aooiaracy  of  determination  of  b  and  o  with  t  .  In  this  example  it  was 
assiimed  that :  -  ® 

(i)  6t  was  0.01  seo, 

(ii)  s  was  of  full  scale. 
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(iii)  the  ncminal  value  cf  b  was  135  secs 
and  (iv)  the  ncminal  value  cf  o  was  15  sec. 

All  cf  these  values  were  ocnsidered  tc  be  reascnably  typical  of  a  medium 
range  ballistic  missile  such  as  Blue  Streak.  Three  separate  cases  were 
considered,  as  follows 

(a)  'a'  unknown, 

(b)  ’a'  known  to  a  standard  deviation  of  IJfi  of  f\all  scale, 


and  (c) 


•a'  known  precisely. 


For  (a)  and  (b)  above,  curves  for  t  -  t  equal  to  10  secs  and  20  secs  are 

'  no 

given.  For  (o),  a  oxirve  for  t^  -  t^  equal  to  20  sees  only  is  given. 


It  will  be  seen  from  the  curves  given  in  these  figures  that  the 

standard  deviations  of  b  and  o  decrease  rapidly  as  t  is  increased. 

0 

Similarly,  for  constant  t^,  the  standard  deviations  decrease  even  more 

rapidly  as  t^  -  t^  is  increased.  These  curves  shew  that  accurate  values 

of  b  and  c  will  only  be  obtained  with  large  values  of  t  and  t  -  t  ,  that 
is  towards  the  end  of  the  boost  phase.  The  analysis  is°also  restricted  to 
the  later  stage  of  the  boost  phase  by  the  assumption  of  zero  aerodynamic 
forces.  This  restriction  means  that  the  analysis  will  only  be  applied  to 
times  when  the  acceleration  is  changing  rapidly  and,  as  a  consequence,  when 
the  mass  flow  and  thrust  are  most  likely  to  be  ohan^ng.  The  effect  of  the 
chan^ng  pump  inlet  pressures  (arising  from  the  changing  missile  accelera¬ 
tion)  on  the  mass  flow  and  thrust  can  be  predicted  using  the  known,  nominal 
engine  characteristics.  The  period  of  the  analysis  must  then  be  chosen  so  that 
the  mass  flow  and  thmist  are  constant  to  within  approximately  C.2^.  This  per¬ 
missible  period  of  analysis  will  be  large  when-  the  missile  has  just  left  the 
atmosphere  and  will  decrease  tov^ards  cut-off.  In  practice,  therefore,  the 
accuracy  is  likely  to  be  reasonably  constat  over  the  later  stage  of  boos'ted 
flight,  the  gain  in  accuracy  due  to  increasing  t^  being  balanced  by  the 

decrease  in  accuracy  due  to  the  reduction  ii^  t^  -  t^. 


If  'a'  is  known  precisely,  a  considerable  improvement  in  accuracy  may 
be  realised  compared  with  that  obtaining  when  'a'  is  linknown.  For  instance, 
in  the  particular  example  given  in  Figs. 8  and  9,  when  t^  =  I50  sees,  and 

t^  -  t^  =  20  sees  a  precise  knowledge  of  'a'  improves  the  accuracy  by  a 

factor  of  8  for  b  and  16  for  o  compared  with  the  case  when  'a'  is  unknown. 

In  practice,  'a*  can  never  be  known  pi^cisely  even  with  very  accurate  pre¬ 
flight  calibration  because  any  bias  errors  arising  in  the  telemetry  system 
dtiring  flight  must  inevitably  be  inoliided  in  the  coefficient  'a'.  If  ’a' 
is  known  to  a  standard  deviation  of  of  full  scale,  a  general  improvement 
in  aocurapy  will  result  compared  with  the  case  of  'a'  unknown  as  is  shown  by 
the  curves  in  Figs. 8  and  9.  The  main  effect  appears  to  be  that  the  sensi¬ 
tivity  of  the  accuracy  to  t^  is  reduced.  It  should  be  noted,  however,  that 

for  the  larger  values  of  t^,  the  two  curves  cross  over  so  that  an  analysis 

assuming  that  ’a’  is  unknown  leads  to  better  accuracy  than  that  when  'a'  is 
known  only  to  I/0  standard  deviation.  This  is,  of  course,  due  to  the  fact 
that  in  these  oircumstanoes  the  method  of  analysis  determines  *a'  to  a 
standard  deviation  of  better  than 
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If  p  is  defined  as  the  inverse  of  the  time  intejrval  between  suooessive 
readings,  6t,  then  p  will  be  equal  to  the  numbers  of  readings  per  seo,  or, 
as  it  is  usually  termed,  the  sampling  frequenoy.  Replaoing  6t  by  p  in 
equations  (25),  (26)  and  (2?)  shows  that  the  standard  deviations,  s^,  s^ 

and  s^  are  proportional  to  the  inverse  of  the  square  root  of  the  sampling 

frequenoy.  Pigs. 8  and  9  were  prepared  assuming  a  sampling  frequenoy  of 
100  readings  per  seo  which  roughly  corresponds  to  the  switch  speed  of  current 
telemetry  systems.  The  maximum  accuracy  will  be  obtained  if  all  of  the 
available  readings  are  used.  However,  as  the  use  of  high  speed  computers  is 
envisaged,  the  numbers  of  readings  employed  in  the  calculation  is  relatively 
vmimpoirtant . 

The  coefficient  o  appears  only  in  the  expression  for  s^,  the  standard 

deviation  of  b.  To  a  limited  extent,  o  (=  Ig/k)  may  be  adjusted  by  the 
choice  of  k,  the  accelerometer  scale  factor.  As  s^  is  inversely  proportional 

to  o,  it  follows  that  k  should  be  chosen  as  small  as  possible  in  order  that 
s^  should  be  as  small  as  possible.  Normally  only  one  accelerometer  is  used, 

its  range  being  chosen  to  cover  the  whole  range  of  accelerations  which  are 
likely  to  ooour  in  the  boost  phase.  This  imposes  a  lower  limit  on  the  choice 
of  k.  The  use  of  two  or  more  aooeleraneters  to  cover  consecutive  phases  of 
the  boost  phase  enables  a  smaller  valtie  of  k  to  be  chosen  but  this  will  only 
have  advantage  if  the  errors  introduced  by  missile  vibration  are  not  iiioreased 
(increasing  s)  as  k  is  reduced.  The  accelerometer  scale  factor,  k,  will  be 
measured  by  the  pre-flight  calibration  of  the  system.  An  estimate  of  the 
standard  deviation  of  k  must  also  be  made.  The  estimated  standard  deviation 
of  the  determination  of  the  specific  impulse  will  therefore  be  the  vectorial 
am  of  the  standard  deviations  of  o  and  k. 

9  CONCLUSIONS 


The  method  as  presented  provides  an  accurate  method  of  determining 
the  vaoum  specific  impulse  of  a  missile  from  longitudinal  acceleration 
records  provided  that  the  engine  performance  is  reasonably  constant.  In 
order  to  achieve  as  high  an  accuracy  as  possible,  6t  should  be  chosen  as 
small  as  possible,  that  is  the  sampling  frequenoy  should  be  as  high  as 
possible,  (t^  -  t^)  the  period  of  flight  to  be  analysed  should  be  diosen 

as  large  as  possible  consistent  with  the  assmption  of  constant  mass  flow 
and  k  the  scale  of  the  accelerometer  should  be  chosen  as  small  as  possible. 
The  accuracy  of  the  method  in  the  case  where  ’a’,  the  zero  bias  of  the 
accelerometer  is  known  is  very  much  better  than  that  obtaining  in  the  case 
where  ’a'  is  unknown,  'k’  should  be  measured  as  accurately  as  possible  by 
pre-flight  ground  testing.  The  method  is  useful  in  that  it  provides  an 
overall  check  on  the  performance  of  the  rocket  engine  and  unlike  other 
methods  does  not  assme  knowledge  of  the  engine  operating  oharacterlstios 
other  than  that  the  engine  is  running  steadily. 

The  method  oan  be  used  directly  to  derive  an  accurate  measurement  of 
the  parameter  m^fi.  Combining  this  with  direct  measttromelit  of  mass  flow, 

e«g.  by  flovT-meter,  oan  give  a  better  measurement  of  missile,  mass  in  the 
later  stages  of  boost-  flight  than  could  be  obtained  by  integration  of  mass 
flow  from  launch. 


Fundamentally,  there  is  no  reason  why  the  mass  flow  or  the  speoifio 
impulse  need  be  assvuned  constant  except  to  simplify  the  analysis.  To  extend 
the  method  to  cover  the  case  when  the  engine  performance  is  varying,  it  is 
only  necessary  to  assume  that  the  mass  flow  and  the  speoifio  impulse  oan  be 
expressed  as  polynomials  in  terms  of  t,  the  time  variable.  The  method  of 
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least  squares  fitting  hy  the  process  of  differential  adjustment,  although 
more  complicated,  will  then  evaluate  the  coefficients  cf  the  polynomials. 
Extension  of  the  method  on  these  lines  would  enable  the  method  to  apply  to 
missiles  such  as  Blue  Streak  which  employ  a  system  of  thrust  reduction 
dijring  the  later  stages  of  the  boost  flight. 
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FIG.  6.  VARIATION  OF  WITH  THE 

LOCAL  MASS  RATIO,  'a’  SUBJECT  TO  ERROR  8a. 


fwfp  3963 


TN.q.W.  S67 

FIG.  7. 
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accuracy  of  determination  of 


T.N.  q.W.  SS7 

FIG.  8 


FIG.  8.  VARIATION  OF  THE  ACCURACY  OF 
DETERMINATION  OF  'b’  WITH  TIME  OF 
START  OF  ANALYSIS. 
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